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Abstract. We discuss the exotic properties of the heat-trace 
asymptotics for a regular-singular operator with general boundary 
conditions at the singular end, as observed by Falomir, Muschietti, 
Pisani and Seeley in [FMPS03] as well as by Kirsten, Loya and 
Park in [KLP06]. We explain how their results alternatively follow 
from the general heat kernel construction by Mooers [Moo99], a 
natural question that has not been addressed yet, as the latter work 
did not elaborate explicitly on the singular structure of the heat 
trace expansion beyond the statement of non-polyhomogeneity of 
the heat kernel. 



1. Introduction 

In this paper we revisit the discussion of the heat kernel and the un- 
usual properties of the heat trace expansion for a general self-adjoint 
realization in L 2 (0, 1) of the regular-singular operator (M + = (0, oo)) 

A, = + 1 (u 2 - i) : C7 ~(R+) CS°(R + ), v > 0, 

where C^°(1R + ) denotes the space of smooth functions with compact 
support in IR + . 

The heat trace expansion of A u for v > is well-understood in 
[KLP08] and in [Moo99], and in fact does not exzibit new phenom- 
ena, compare in particular the general discussion of Gil, Krainer and 
Mendoza [GKM10]. The intricate case is rather v = 0, where the rela- 
tion between the explicit approaches of [FMPS03], [KLP06] and the 
methods of the heat kernel construction in [Moo99] is less obvious. 

The regular-singular operators A„ arise naturally in the spectral geom- 
etry of spaces with isolated conical singularities, modelled by a bounded 
generalized cone (^(iV) = R + x N,g = dx 2 © x 2 qn) over a closed Rie- 
mannian manifold (N, g ) of dimension n. This includes the example 
of a higher-dimensional disc 3) C M ra+1 , where the euclidean metric 
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takes the form dx 2 © x 2 gs^ with respect to standard polar coordinates. 
The Laplace-Beltrami operator on (^(N) = M + x N, g) is a symmetric 
operator in L 2 (^(iV), g) and takes the form 

Ay(jv) = -9,. - - <9 X + X _2 Aiv. 
Under the unitary transformation (see [BrSe87]) 

$ : L 2 ($L + ,L 2 (N,g N ),dx) ->■ L 2 (tf (N) , g) , u ^ x n/2 co, 
the Laplacian A<^tm unitarily transforms to 

$ o A w o = -«9 2 + 1 ^Ajv - =: Aj (iV) . 

The spectral decomposition on the cross section (N,g^) decomposes 
A^jyA into a direct sum of operators A„ over the z/ 2 -eigenspaces of A N . 
We are interested here in the operator over the kernel of Ajv 

d 2 

dx 2 Ax 2 

The spectral properties of A for general self-adjoint boundary con- 
ditions have been subject of careful analysis by Falomir, Muschietti, 
Pisani and Seeley in [FMPS03] as well as by Kirsten, Loya and Park 
in [KLP06]. In both instances the authors uncovered new unsual phe- 
nomena in the heat trace expansion of A, if the boundary conditions 
at x = do not define the Friedrichs extension. Their (and our) result 
reads as follows. 

Theorem 1.1. The self-adjoint boundary conditions of A at x = are 

parametrized by 9 G [0, n), where 9 = tt/2 corresponds to the Friedrichs 
extension. We write A(8) for the corresponding self-adjoint realization. 
With 

'■= 1 — log 2 + tan 9, 

where 7 denotes the Euler-Mascheroni constant, we find for 9 7^ 7r/2 



A = -—- — :C™(R+)^C~(I 



„i 00 
Tr ie - tA W:= / e- tA W(x,x)dx~T Tl e- tA ^ + V % t' 
Jo j=0 



+ -Im / (logy + in + 2n e )~ 1 dy , t -> 0. 

7T \Ji y J 

Both discussions of [FMPS03] and [KLP06] have been performed in- 
dependently from the earlier work by Mooers [Moo99], who in partic- 
ular constructed the heat kernel of A for a general self-adjoint realiza- 
tion and observed its non-polyhomogeneity in case of a non-Friedrichs 
extension at x — 0. 
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However, [Moo99] did not elaborate further on the particular structure 
of the non-polyhomogeneous heat kernel and its heat trace behaviour, 
leading to the natural question of whether her analysis can be reconciled 
with the explicit results in [FMPS03] and [KLP06]. Hereby, certain 
inconsistencies in [Moo99] become apparent, see Remark 4.2. 

In this note we address this issue and show how the heat trace expan- 
sion results in [KLP06] follow straightforwardly from the heat kernel 
construction of Mooers in [Moo99]. At various points we chose to 
provide arguments, alternative to [Moo99]. 

This paper is organized as follows. We first classify the self-adjoint re- 
alizations for A in §2, and study the solution to the signaling problem 
in §3. The solution to the signaling problem is the central ingredient 
in the construction of the heat kernel for general self-adjoint bound- 
ary conditions, which is explained in §4 and is basically a revision of 
[Moo99]. In §5 we establish well-definement of the various expres- 
sions in the general heat kernel formula. In the final §6 we derive the 
heat trace expansion directly from the heat kernel structure in §4, thus 
reproving Theorem 1.1 by a different method. 



The author would like to thank Thomas Krainer for helpful remarks 
and gratefully acknowledges the support by the Hausdorff Research 
Institute at the University of Bonn. 
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2. Self-adjoint boundary conditions 



Consider the differential operator 
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acting on C£°(R + ) C L 2 (R + ). We always denote by Cg°(I) the space 
of smooth functions with compact support in I C R. We define the 
minimal closed extension A min to be the graph closure of A, which 
is a densely defined and symmetric. The maximal closed extension is 
A max . (A mm ) . 

We then have the following the well-known characterization of the 
maximal domain 5?(A max ), written out in various sources, including 
[BrSe87], [Moo99], [FMPS03],[KLP06], [LeVeII] and [Ver09]. 

Proposition 2.1. Let f G ^(A max ). Then there exist constants 
c + (f), c_(/) G C depending only on f , and a continuously differentiable 
/ G ^(A min ) ; with J{x) = 0(x 3 / 2 \og(x)) and f'(x) = 0{x x l 2 log(x)), 
as x y ; such that 

f{x) = C+(f)y/x + C-(f)y/xl0g(x) + f(x). 

Moreover, for any f,g G ^(A max ) the following Green's identity holds 



(Af,g) - (f,Ag) = c4f)c + (g) - c + {f)c„{g) 



An extension i^(A mm ) C @ C ^(A max ) is self-adjoint, if 

9 = {/ G ^(A max ) | V 5 G ^ : (Af,g) L2 = (f,Ag) L2 }. 

Consequently we obtain a self-adjoint extension of A by choosing the 
following boundary operators on / G i^(A max ) 

B e (f) := cos9 -c + (f)+ sm0 -c-(f), 8 G [0, 7r). (2.1) 

We define then 

9{A{8)) := {/ G ^(A max ) | £ fl (/) = 0}. (2.2) 

The extensions A(8),8 G [0, tt) define self-adjoint realizations of A by 
Proposition 2.1, and in fact classify completely all self-adjoint boundary 
conditions at x = 0. 

3. The signaling solution 

The fundamental component in the heat kernel construction of Mooers 
in [Moo99] is the signaling solution F(h)(-,t) G @(A max ),t G Rq, 
with R^ = [0, oo), defined for any given h G Aoc O^o ) n C (R + ), as a 
solution to the so-called signaling problem 
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{d t + A)F(h)(x,t) = Q, F(h)(;0) = 0, 

C„(F(h)(;t)) = h(t), t>0. [ ' 

Note that there does not exist 9 G [0, tt), such that F(h)(-, t) G @(A(6)) 
for t G Kq, since by uniqueness of solutions to the heat equation, 
F(h)(-, 0) = then implies F(h) = 0. The signaling solution is in fact 
unique, since for c-(F(ti)(-,t)) = 0, F(h)(-,t) G @(A(n/2)) and hence 
F(h)(-, 0) = then implies F(h) = 0. 

Let the differential expression A act on C^°(1R + ) and consider the heat 
kernel E + (x,x,t) of the Friedrichs extension of A in J 2 (R + ). We put 
N'E + := x~ 1/2 E + (x } x } t) \ x = 0. By [Les97] and the asymptotics of 
the modified Bessel functions, see [AbSt92] 



I viz) ~ ; r, as Z — > 0, 



we have 



. _ o . — -Q 

/yrp I nn ff \ / ^ I ^ 



E + (x,x,t) = — — J — exp 



2t \2t J V 4t 



(3.2) 



Lemma 3.1. [Moo99, Proposition 4.2] 

For any h G J 1 1 oc (IRq ) fl C°°(1R + ) the signaling solution is given by 

F(h)(x,t) :=- f h(t-T)N'E+(x,r)dT. 
Jo 

Proof. [Moo99, Proposition 4.2] establishes the statement by deriving 
the formula for F(h) conceptually. Here, for purposes of brevity, we 
simply show by a direct computation that the solution F(h) above 
indeed provides the signaling solution. With N'E + solving the heat 
equation, clearly (d t + A)F(h)(x,t) = 0. 

By the explicit formulas (3.2) we find for any fixed t > 0, F(h)(-,t) G 
J 2 (M + ). Similarly, differentiating F(h)(x,t) explicitly in the first ar- 
gument, we find AF(h)(-,t) G J 2 (R + ) and consequently, F(h)(-,t) G 
^(A max ). 

It remains to identify c-(F(h)(-,t)). We write h(t-r) =: h(t)-rJF(t- 
r), where Jtf G C°°(1R + ) by smoothness of h. Then, substituting T = 
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x 2 jr in the integrals below, we find 



F(h)(x,t) = -[ h(t-r)N'E + {x,T)dT 
Jo 




=: h + h- 

For the first integral and x 2 < t (note t > is fixed) we obtain 




= h(t)y/xlog(x) + 0(y/x), X -»- 0. 

For the second integral I2 we obtain similarly I2 = 0{yfx) as x — > 0, 
and hence indeed find c-(F(h)(-,i)) = h(t). □ 



It should be noted that the naive expansion of F(h) does not lead to a 
closed expression of c + (F(h)), which requires rather integral transform 
arguments and is the fundamental aspect of the next section. 

4. Heat kernel for general boundary conditions 
We can now discuss the heat kernel for a self-adjoint operator 
0(A(0)) = {/ G @(A max )\B e (f) := cos 6 c + (f) + sin£c_(/) = 0}, 

with 9 G [0,7r). By the Duhamel principle, the heat trace expansion 
of any self-adjoint extension of A in L 2 (0, 1) with boundary conditions 
Bq at x = and separated boundary conditions at x = 1, differs from 
the heat trace expansion of A(9) (trace integral taken over [0,1]) by 
the classical trace contribution from the regular boundary x = 1 and 
0(t°°),t -> 0, cf. for instance Lesch [Les97, Theorem 1.4.11]. 

The heat kernel E + for the Friedrichs extension A(n/2) is well- 
known with the explicit expression written out in (3.2). Fix any 
6 G [0, 7r)\{7r/2}. We do not distinguish notationally between the 
pseudodifferential operators and their Schwartz kernels. Consider 
4> G Qj°(Ro )> vanishing to infinite order as x — > 0, and put u = E + (f>. 
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We seek to correct u by F{h) for an appropriate h G L 1 1 oc (Mq )nC°°(IR + ), 
to satisfy the boundary conditions B e {u + F{h)) = 0. Note 

uj :=u + F(h) G 0(A(0)) 
& c+(u) + c+(F(h)) = -tan(0)c_(F(/O) (4.1) 
^ + G * t h(t) = - taa(0) h(t), 

where G is the convolution kernel mapping h = C-(F(h)) to c + (F(h)). 
We explify G(£) using the Laplace transform. For any g G ^i oc (K^") D 
C 0O (M + ), not growing exponentially as £ — >■ oo, the Laplace transform 
is defined as follows 

{S?g){()= [ ^(£)exp(-C£)rf£,Re(C)>0. (4.2) 

The inverse Laplace transform is given for any any 5 > and analytic 
L((), integrable over Re(£) = S, by 

{J?- 1 L){t) = ±, f e«L(()d{. (4.3) 

Lemma 4.1. 

^fC(C) = logVC + 7-log2. (4.4) 

Proof. We compute 

(SfF(h))(x,0 = -(-270(C) ■ (J?N'E + )(xX) 

= -Vx~K {x^()(J?h)((), Re(C)>0 

where Kq is the modified Bessel function of second kind, and in the 
definition of a/C we fi x the branch of logarithm in C\M _ . We assume 
that h(t) is not of exponential growth as £ — > oo, so that 
is well-defined for Re(£) > 0. The Bessel function K (z) admits an 
asymptotic expansion, see [AbSt92] 

K (z)~-log(z) + (log2- 1 )+K (z), K (z) = O(z), z^O, (4.5) 
where 7 6 R is the Euler constant. Consequently 

c + ((J?F(h))(; 0) = -(^)(C)(log2 - log VC - 7), f4 6) 

c.((J2fF(/i))(.,C)) = (^)(C). 
Taking the inverse Laplace transform, we obtain 

F(h)(x,t) = v / ilog(x)^f- 1 (c_(^fF(/i))) + v^ _1 (c+(.£fF(/i))) 

+ y/ESf-\(Sfh)K (xyP)), 
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where each Se- x (c±(J?F(h))) exists and ^{{^KjK^xyf)) = 
0(x 3//2 ), uniformly as x — » 0. Consequently, indeed 

S?(c±(F(h))) = c±(S?F(h)). 

This yields the following explicit expression for the Laplace transform 
ofG 

(^G)(C) = lo gv ^+7-log2. (4.7) 

□ 

As already noted at the end of §3, the explicit expression for c+(F(h)) = 
G *t h(t), obtained in Lemma 4.1 by means of integral transforms, 
cannot be obtained by naive expansion of F(h)(x,t) in x directly. 

Remark 4.2. Similar argument may be performed using Fourier trans- 
formation instead of the Laplace transform. The Fourier transform has 
been used in [Moo99, Lemma 4.4], where however the factor (7 — log 2) 
from (4.7) is incorrectly missing. More precisely, the scaling argument 
outlined in [Moo99] determines the Fourier and the Laplace trans- 
form of G to be given by log y/^ only up to an additive constant, which 
cannot be specified by her method. 

Returning back to (4.1) and taking the Laplace transform of its third 
identity, we obtain 

h{t) = -if" 1 {(££G) + tan^))- 1 *t N'E + <P =: - Kg * t 

u + F{h) = (E + + N'E + * t Kg * t N'E+) e 0(A(0)), ' 

where the inverse Laplace transform in the formula for Kg is given 
explicitly as follows. For any 5 > and Kg 1=7 — log 2 + tan 9 we have 
a priori 

Kg(t) = / -. -= ^ d( = / ~, -= r- d( 

2« J s+iR (log VC + Ke) 2m JiR (log VC + Kg) 



- lim — / e itR e xt flog y/R + iarctan(x/5)/2 + Kg) dx 
R-^oo 2ni J V / 

+ lim — / e~ itR e xt flog VR- z arctan(x/5)/2 + Kg) dx. 
R-^-oo 2%i Jq V / 

The latter two integrals behave as 0((log-R) _1 ), as R — > 00, and hence 
vanish in the limit. Consequently 

Kg(t) = l^-f m e* C (log + 7 - log 2 + tan(0)) _1 d(. (4.9) 
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Well-definement of Kg* t N'E + <j), and regularity of h, F(h) is not obvious 
and requires a detailed analysis of Kg, which is the content of §5 below. 
Then, in view of (4.8) we have proved the following 

Theorem 4.3. The heat kernel Eg ofA(9),6 G (0, it) is given by 
E e = E + + N'E + * t Kg * t N'E+. 



5. Analysis of the convolution kernel Kg 
The first key step here is a specific integral representation of Kg(t). 



Proposition 5.1. Let k{6) := 7 — log 2 + tan(#). Then there exists a 



bounded Kl G C°°(E£), such that 



K e(t) = |^°° e-*» (log tfj + in/2 + k{6))- 1 | + K l e {t) 

/oo 
e-^ ((logy + 2k{6)) 2 + n 2 )- 1 dy + K}{t). 

Proof. We compute from (4.9) 

Kgit) = 5^ e* (lo gv ^+^)) _1 d( 

= ±-( e^Uog^C + ^O)) 1 d( 

+ — / e itx (log + W4 + _1 dx 

1 /"°° / \ _1 

+ — J e~ ltx (log y/C-iir/4 + k(9)) dx 

=:Kl{t)+K 2 e {t) + Kl(t). 



The first summand Kg G C°°(]Rq~) is smooth and bounded. For the 
summand Kg we deform the integration contour to i[l,oo). Let the 
contour 7^ := {Rexp(i(f>) \ (f> G [0, 7r/2]} be oriented counterclockwise. 
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R 



FIGURE 1. The integration contour 7^. 



We use the O-notation for the asymptotics as R — > 00. Then, substi- 
tuting x = Rexp(i(p), we find for some C > 



/ 



e itx (log y/x + in/A + k(9)) 1 | dx 

n f n/2 . exp (-tR sin + .. , 

= R I — ; 7^ 777; — j^— I Ct0 

./() 



- R l 

Jo 



log VR + i((f>/2 + Tr/4) + K(9) 
w/4 exp(-tR sin <f) + i(<j) + tR cos <f))) 



d(f) 



log Vi? + i(0/2 + tt/4) + «(0) 

<- / exp(-t J Rsin0)cos0# + O( J R"°°) 

log/? Jo ' 

where in the fourth line we have used the fact that cos is bounded from 
below for e [0, 7r/4], and in the final line we substituted y = sin0. 
Consequently, we may indeed deform the integration contour of Kg to 
00) and find 

• poo 

K 9<t) = to / e"* (tog Vy + <7r/2 + /6(e))" 1 dy 
+ / e ite (log v / a : +^/4 + k(6'))~ 1 dx, 

where the second summand is clearly smooth and bounded on Rq , 
since Im(x) > for x e 71. We denote the second summand in the 
expression for Kg by Kg e C°°(Rq). By a similar exercise we deform 
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the integration contour of Kg to (— oo)) and find 

■ poo 

K$(t) = (log y/y - m/2 + k(9))- 1 dy 

+ f e~ itx (\og^ - in / A + k{6)Y 1 dx, 

where 7! = {z \ z G 71} is oriented clockwise; again the second sum- 
mand is clearly smooth and bounded on Rq", since Im(x) < for 
x G 7i- We denote the second summand in the expression for Kg 
by Kg G C°°(Ro ). In total we now obtain using i(z — z) = —2lm(z) 
and setting K\ := K] + K$ + K\ 

K e(t) = {^°° e~ ty (log + ztt/2 + k^))" 1 } + 

/oo 
e"* ((logy + 2/<£)) 2 + vr 2 )- 1 dy + tf*(t). 

(5.1) 
□ 

As a consequence of Proposition 5.1, we find that Kg G Li oc (Rq") fl 
C°°(R + ) is bounded as t — > 00. Integrability at t — follows using 
Fubini theorem. Indeed we find for any T > 0, using (5.1) 

/oo pT 
J e~ t y((\ogy + 2K(6)) 2 + n 2 )- 1 dy 

r°° f\ _ e - Ty \ 1 

= ^ f J ((log?/ + 2 K ^)) 2 + 7T 2 )- 1 dy < OO. 

On the other hand, also N' E + <p G Aoc^o ) n <^°°(R + ), bounded as 
t — )• 00. Indeed, substituting X = x/y/t we find for supp(</>) C [0, 1) 

N'E +( f>(t) = J q 1} § exp (-^Q <f>(x) dx 

= r 1 " exp ftXVt) dX, 

which shows integrability at t = 0. Consequently, the convolution 
h = Kg* t N'E + 4> indeed exists and h G C°°(R + ) is continuous at t = 
and not growing exponentially as t — > 00. In particular, the Laplace 
transform J^(h)(() is well-defined for Re(£) > 0. Similarly, for each 
fixed x G Rq , the signaling solution F(h)(x, •) G C°°(R + ) is continuous 
at t = and its Laplace transform JzfF(/i)(£) is also well-defined for 
Re(C) > 0. 
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6. Unsual heat trace expansions 

In this final section we study the asymptotic expansion of Tri(Eg), 9 G 
[0, 7r)\{7r/2} and establish Theorem 1.1. Discussion of the kernel 
N'E + * t Kq * t N'E + is central and consists of two steps. We begin with 
an explicit evaluation of asymptotics of TriN'E + * t N'E + as t — > 0. 

Proposition 6.1. 

r 1 i 

Tri N'E + * t N'E + = / (N'E+ * t N'E+)(x, t) dx = - + 0(t°°), t -> 0. 
Jo 2 

Proof. We first substitute y = x 2 

[ {N'E + * t N'E + ){x,t)dx = [ [ N'E + {x,t-r)N'E + ){x,T)dTdx 
Jo Jo Jo 



1 rt X ( 2 / 1 1 



, , . exp —x — — — dr dx 

, J 4r(t-r) P V \4(t-r) 4r ' 1 

exp ( -y ( — — — r - — ] ) drdy 



JO 



(t - r) " V V4(* " t) 4r 



By Fubini theorem we may interchange the integration orders and find 
after integration in y 



1 '* 1 / / 1 1 



ex P _ 2/ tt: 7 -7- rfr ^ 



o 7 8r(t-r) " V \4(t - t) At 

1 y* / / * \\ . 1 



where in the last step we have estimated for r G [0, t] 
6XP (-47^) " GXP ("41 



We can now discuss the trace of full kernel N'E + * t Kg * t N'Ej, 
Theorem 6.2. 



□ 



Tri (iV'E+ * t K e * t N'E+) ~ 6^' 



i=o 



1 / Z" 00 e~ tJ/ \ 
+ -Im(/ — -(log v ^ + i7r + 2K(^))J dy, 
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Proof. Using Proposition 5.1 we find 
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Tn (N'E+ * t K e * t N'E+) 

= 2 {N'E + * t N'E + ){x,r)dydrdx 

Jo Jo J i {{logy + 2k{6)) 2 + it 2 ) 

+ I [ K]{t — r){N'E + * t N'E + ){x, r) drdx =: T\ + T 2 . 
Jo Jo 

We want to apply Fubini theorem and interchange the integration or- 
ders. For this we note using Proposition 6.1 

OO ft pi p -{t- T )y 

{N'E + * t N'E + ){x, t) dx dr dy 



1 Jo Jo 



{{\ogy + 2K{6)) 2 + 7t 2 ) 

(\ \ f 00 e~ ty -l 

due to log(y) 2 behaviour in the denominator of the integrand. Moreover 

I [ K]{t-T){N'E + * t N'E + ){x 1 T)dxdT 
Jo Jo 

= jfjfJ(f-T) Q + O(o) dr 

< Q + °( t0 °)) f^Kl{t-r)dr<^ 

due to smoothness of Kg. Consequently we may indeed interchange the 
integration orders and obtain using Proposition 6.1 and the integral 
representation of Kg{t) in Proposition 5.1 

^ = jH J Im { e -<t- T » (log y/y + in/2 + ^fl))" 1 } 

x f {N'E + * t N'E + ){x,r)dxdTdy (6.1) 
Jo 

11 i m ) ^ tl — y dy + (f 

^ y{\ogjy + m + 2K{6)) ( 
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as t — > 0. For T 2 we obtain 



T 2 = f K](t - r) [ {N'E + * t N'E+)(x, r) dx dr 
Jo Jo 



(6.2) 



since ifg G C°°(Rq"). The expansions (6.1) and (6.2) together yield a 
full asymptotic expansion of Tri (N'E + * t Kq * t N ; E + ) 

Tn (N'E + *, K, *, N'E + ) 



iut— 

7T \Ji y 



{log^/y + iir + 2 K (9)) dy 



3=0 



□ 



Together with Theorem 4.3 we obtain the statement of Theorem 1.1. 
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